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ABSTRACT 

We obtained Hubble Space Telescope /Wide Field Camera 3 imaging of a sample of ten of the nearest 
and brightest nuclear clusters residing in late-type spiral galaxies, in seven bands that span the near¬ 
ultraviolet to the near-infrared. Structural properties of the clusters were measured by fitting two- 
dimensional surface brightness profiles to the images using GALFIT. The clusters exhibit a wide range 
of structural properties, with F814W absolute magnitudes that range from —11.2 mag to —15.1 mag 
and F814W effective radii that range from 1.4 to 8.3 pc. For six of the ten clusters in our sample, we 
find changes in the effective radius with wavelength, suggesting radially varying stellar populations. In 
four of the objects, the effective radius increases with wavelength, indicating the presence of a younger 
population which is more concentrated than the bulk of the stars in the cluster. However, we find a 
general decrease in effective radius with wavelength in two of the objects in our sample, which may 
indicate extended, circumnuclear star formation. We also find a general trend of increasing roundness 
of the clusters at longer wavelengths, as well as a correlation between the axis ratios of the NCs and 
their host galaxies. These observations indicate that blue disks aligned with the host galaxy plane 
are a common feature of nuclear clusters in late-type galaxies, but are difficult to detect in galaxies 
that are close to face-on. In color-color diagrams spanning the near-UV through the near-IR, most of 
the clusters he far from single-burst evolutionary tracks, showing evidence for multi-age populations. 

Most of the clusters have integrated colors consistent with a mix of an old population (> 1 Gyr) and a 
young population 100-300 Myr). The wide wavelength coverage of our data provides a sensitivity 
to populations with a mix of ages that would not be possible to achieve with imaging in optical bands 
only. The surface brightness profiles presented in this work will be used for future stellar population 
modeling and dynamical studies of the clusters. 

1. INTRODUCTION 

Hubble Space Telescope (HST) imaging surveys have 
shown that most galaxies contain bright, compact stel¬ 
lar systems known as nuclear clus ters (NCs) in their 
photometric and dynamical centers (iCarollo et al.l 119971: 
iMatthews et all 119991: iBoker et~ahl l2QQ2h . Similar in 
many respects to high-mass globular clusters (GCs), NCs 
typically have velocity dispersions of 15 - 35 km s _1 , ef¬ 
fective radii (or half-light radii) of a few parsecs , and 
dynamical masses of 10 6 - 10 7 M 0 (Bo ker et al.l H99l 
l2002t IWalcher et ~aJl l200Rh . Unlike GCs, NCs in spiral 
galaxies have complex star formation histories, undergo- 
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ing repeated bursts of star formation, with spectra of- 
ten consistent with a continuous star formation histor y 

(|Walcher et al. 1(20061: iRossa et~ahll2QQ6l: ISeth et akll2QQ6h . 

Young (<100 Myr) populations are found in most spiral 
NCs, and there is limited evidence for complex p opula- 


tions in early-type nuclei as well ( 

iMonaco et al.l 2009: 

Seth et al. 20101: Lvubenova et alj 

2013 

). In a recent 

HST archival studv. iGeoredev & Boker ( 

20141) analvzed 


the photometric and structural properties of a sample of 
228 NCs in spiral galaxies, and confirmed that most NCs 
are similar in size to GCs, but found that the largest 
and brightest NCs in spiral galaxies occupy the regime 
between ultra-compact dwarfs (UCDs) and the nuclei of 
early-type galaxies in the size-luminosity plane, which is 
consistent with the idea that some UCDs and GCs are 
the remnant nuclei of tidally disrupted spiral g alaxies 
(|Bokerll2008l : iMieske et al.l 120131: ISeth et ahll20141 ). They 
also found that the effective radius of NCs tends to be 
smaller when measured in bluer filters, suggesting that 
younger stellar populations tend to be more concentrated 
than the rest of the stars in NCs. We note that the 
Milky Way’s NC also has a popu lation of young stars 
concentrated in the central 0.5 pc (jPaumard et al1 l2006: 
iLu et al.l 120131: iSchodel et al.ll2014] h 

Two main scenarios have been proposed to explain 
the formation of NCs: in-situ formation from the ac¬ 
cretion of l ow angular m omentum gas onto the galac¬ 
tic nucleus (|Milosavlievicll2004} ), and mergers of globular 
clusters migrating t owards the galac ti c nucleus due to 
dynamical friction (iLotz et al.l l20(Vfl L lHartmann et all 
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( 2011 !) compared simulations of the migratory/merger 
scenario with integral field unit (IFU) observations of 
NGC 4244, and found that in-situ formation accounts for 
at least 50% of t h e NC’s mass. lAntoninil ( 2013L120141) and 

lArca-Sedda fe Capuzzo-Dolcettal (120141 ) demonstrated 
that the migratory/merger scenario predicts correlations 
between the masses of NCs and global properties of their 
host galaxies that match observed correlations. The 
young populations observed in most NCs rule out the 
sinking of globular clusters as the sole process for NC 
formation (i.e., the sinking of younger star clusters or 
gas accretion must play some role). However, our under¬ 
standing of the formation of NCs is incomplete and the 
reason for the ubiquity of NCs still remains a mystery. 

NCs and BHs are known to coexist in gal axies of all 
Hubbl e types across a wide range of masses (jSeth et al.l 
l2008al) . but some NCs have no central BH down to highly 
constraining limits. The NC in NGC 493 5, for exam¬ 
ple, contains a NC ([Filippenk o fe Holl2003[ ) as well as a 
Type 1 active galactic nuc leus (AG N) powered by ac¬ 
cretion onto a central BH (jFilippenko fe Sargent! 11989 ). 
with mass estimates ranging from 5 x I Q 4 to 3.6 x 10 5 M© 
(jPeterson et al.l[20051: lEdri et al.ll2012[ ). while dynamical 
studies of the NC in M33, a galaxy with a similar mass 
to NGC 4395, have placed remarkably tight upper lim¬ 
i ts of 1500 - 3000 M© on the mass of an y central BH 
(jGebhardt et al.ll200lHMerritt et al.ll2QQlf ). These stud¬ 
ies indicate that at least some NCs can host a central BH, 
but the overall occupation fraction of BHs in NCs is un¬ 
known. The masses of both NCs (Ferrare se et al.l 1 2006: 
Wehner fe Ha rris 200 6|) and BHs dFe rrarese fe Merrittl 
20001: IGehharrlt, et, a.l]l2000l IHaring Rr. Riyll2004ti are cor¬ 
related with the masses of their host galaxy bulges. 
Galaxies with mass above 10 10 M© tend to host BHs, 
while less massive galaxies tend to host NCs, with a tran¬ 
sition region for galax ies with mass 10 8 - 10 10 M© where 
both objects coexist (jGraham fe Spitlerll2QQ9h . Interest¬ 
ingly, for galaxies known to contain both a NC and a cen¬ 
tral BH, the ratio (Mbh + MNc)/Mb u i g e (where Mb u i ge 
denotes the mass of the bulge or pseudobulge compo¬ 
nent of the host galaxy) shows less scatter than eithe r 
Mbh /Mbuige or M N c/M b uige flKormendv fe Hoi 120131) . 
suggesting a link between the building of NCs and BHs. 
These studies add to a growing body of evidence that 
NCs and BHs are both generic by-products of galaxy 
formation, and the gr owth mechanisms of N Cs and BHs 
are somehow related (jFerrarese et a l. 2006). 

Dynamical studies of NCs in late-type spirals have 
been proven useful as a method for stud ying BHs. Pre¬ 
vious studies of I C 342 (Bo ker et al.l 1 1 9991 ) and NGC 
3621 (jBarth et al.l l2009h have constrained A/rh by us¬ 
ing HST imaging and integrated stellar velocity disper¬ 
sion measurements of the NC in order to carry out Jeans 
m odeling of the clusters. Kormendv et al.l (2 0101 ) and 
iNeumaver fe Walcherl (j2012l ) have also performed Jeans 
modeling of NCs in order to set upper limits the masses 
of central BHs. Adaptive optics (AO) assisted IFU ob¬ 
servations can resolve the kinematics of the nearest and 
brightest NCs and increase the sensitiv ity to detect ions of 
BHs within the clusters. For example, iSeth et ah (201 0) 
placed an upper limit of ^ 10 5 M© for the central BH in 
NGC 404 from dynamical modeling of the stellar kine¬ 
matics of the NC using a combination of AO-assisted 


near-IR IFU spectroscopy, optica l spectroscopy, and HST 
imaging. Using similar methods. iHartmann et all ( 2011 ) 
measured a mass of 1.1 x 10 7 M© for the NC in NGC 
4244, and found that if the NC hosts a central BH at all, 
its mass is less than 1% of the NC mass. 

Studies of NCs have mainly been limited by the angular 
resolution of existing telescopes. The next generation of 
large ground-based telescopes, such as the Thirty Meter 
Telescope (TMT), will be able to resolve the kinematic 
structure of stars within the sphere of influence of BHs 
with much lower mass than is currently possible (down 
to ^ 10 4 M©), and will be crucial for future dynami¬ 
cal searches for IMBH within NCs (iDo et al.l 1120141 ). In 
a simulation of observations of NCs with t he E urop ean 
Extremely Large Telescope (E-ELT), iGullieuszik et al.l 
(j2014l) found that observations of NCs with telescopes 
that can resolve individual stars in the cluster will lead to 
great improvements in stellar population and age-dating 
studies, which are currently limited by the various degen¬ 
eracies that come with fitting stellar population synthesis 
models to the integrated cluster light. The nearest and 
brightest NCs will be prime targets for the next genera¬ 
tion of large ground-based telescopes. 

We present a set of uniform, high-signal-to-noise 
( S/N ), unsaturated images of NCs with wavelength cov¬ 
erage from the near-ultraviolet (UV) to the near-infrared 
(IR) for a sample of ten of the nearest late-type galaxies 
hosting NCs. Our data represent a substantial improve¬ 
ment in spatial resolution, wavelength coverage, and ho¬ 
mogeneity over previous imaging studies of NCs. Al¬ 
though all of the galaxies in our sample have archival 
HST images available from previous studies, the archival 
data are very heterogeneous, including images that are 
saturated in the core of the NC and images taken with 
several different cameras. The homogeneity of our data 
allows for a more detailed analysis of the structure and 
stellar populations of NCs than would be possible with 
archival data alone. Detailed structural information of 
the best resolved NCs across a wide range of wavelengths 
may provide clues about formation mechanisms for NCs, 
which are currently poorly understood. 

Here, we describe the detailed two-dimensional (2D) 
structure of this sample of NCs in a uniform way across 
UV, optical, and IR bands. Surface brightness profiles, 
effective radii, magnitudes, and colors are presented. In 
future work, we will describe the stellar populations and 
model the stellar dynamics of the NCs in order to con¬ 
strain the stellar mass and the mass of any central BH 
in each NC. Surface brightness profiles and stellar mass- 
to-light profiles of the NCs will be crucial inputs for dy¬ 
namical models of the clusters. 

The remainder of this paper is organized as follows. 
Section 2 describes the NC sample, the HST observa¬ 
tions and data reduction. In Section 3, we describe sur¬ 
face brightness profile fits to the images of the clusters. 
General results from these fits, which include structural 
parameters and magnitudes of the NCs, trends of cluster 
size and roundness with wavelength, and a comparison 
of simple stellar population models to measured NC col¬ 
ors are given in Section 4. A brief discussion of each 
individual cluster is given in Section 5. Finally, Section 
6 summarizes our most important results and describes 
how they will applied in future work. 
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2. DATA 

We obtained HST/ WFC3 images of ten of the near¬ 
est and brightest nuclear clusters in seven different filter 
bands, spanning the near-UV to near-IR in wavelength. 
The images were taken over a total of 20 HST orbits for 
the GO-12163 observing program. 

2.1. Sample Selection 

The galaxies in our sample were selected according to 
the following criteria: Hubble type in the range Sc-Sm, 
distance < 5 Mpc (with one exception), and the pres¬ 
ence of a nuclear star cluster with magnitude brighter 
than V = 19 mag. The Hubble type selection criterion 
ensures that the galaxies in our sample have little or no 
bulge component. NGC 3621, at a distance of 7 Mpc, 
was included in the sample because it is a rare exam¬ 
ple of a nearby galaxy containing a Type 2 AGN in it s 
nuclear cluster ([Satvaoal et al.l[20071: iBarth et al.ll2009h . 
Basic properties of the galaxy sample are presented in 
Table [U The morphological types, R-band magnitudes 
and heliocentric radial velocities listed in this table are 
from The Third Reference C atalogue of Bright Galaxies 
(jde Vaucouleurs et al.l Il99lh and were looked up in the 
NASA/IPAC Extragalactic Database fl For eight of the 
galaxies in our sample, the distances quoted were mea¬ 
sured using Cepheid variables. We used tip of the red 
giant branch (TRGB) distances for NGC 2976 and NGC 
4244, because Cepheid variable distances were not avail¬ 
able. 

2.2. Description of Observations 

Each cluster was observed in seven filter bands over 
a total of two orbits. For each galaxy, a series of four 
exposures was taken in each filter band. These four ex¬ 
posures were offset by fractional-pixel shifts relative to 
each other using a four-point box dither pattern. The 
final images were constructed by combining these expo¬ 
sures, which allowed for the removal of cosmic-ray hits 
and bad pixels. 

In the UVIS channel, we used the F275W, F336W, 
F438W, F547M, and F814W filters. We used F547M in¬ 
stead of a wide-band filter for the V band because it 
does not include flux from the [O III], Ha, or H /3 emis¬ 
sion lines. The UVIS channel CCD has a pixel scale of 
0704 pixel -1 . At a distance of 5 Mpc, a pixel on the 
UVIS CCD subtends 0.96 pc. The point spread func¬ 
tions (PSFs) range in full width at half maximum from 
07092 in F275W to 07196 in F153M. These angular sizes 
correspond to physical sizes of 2.20 pc and 4.70 pc at 
a distance of 5 Mpc. Since NCs typically have effective 
radii of a few pc, WFC3 provides sufficient spatial reso¬ 
lution for studies of their structure. In order to observe 
each NC with several filters in just two orbits, it is nec¬ 
essary to reduce readout time and eliminate buffer dump 
overheads. To accomplish this, we used the MS512C sub¬ 
array readout mode, which gives a 20" x 20" field of view. 
This field of view is more than wide enough to cover the 
entire NC and sample the region where the NC surface 
brightness profile merges with that of the host galaxy. 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics and Space 
Administration. 


In the IR channel, we used two medium-band filters, 
F127M and F153M, to sample the J and H bands be¬ 
cause the NCs are bright enough to saturate with very 
short exposure times (< 10 seconds) using wide-band 
filters. The IR channel HgCdTe array has coarser res¬ 
olution than the UVIS channel CCD, with a scale of 
0713 pixel -1 . The IRSUB256 subarray was used, giving 
a 35" x 31" field of view for the IR images. We used 
SPARS 10 sampling sequences with NSAMP=10, with 
the exception of the two brightest NCs, IC 342 and M33, 
where RAPID readout mode was used to avoid satura¬ 
tion. 

Throughput curves of the HST/WFC3 filters are 
shown in Figure [lj Table [2] lists the central wavelengths 
and band widths of the filters as well as the exposure 
times used. Shorter exposure times were used for the 
two galaxies with the brightest NCs, IC 342 and M33, in 
order to avoid saturation. Due to the shorter exposure 
times and RAPID sampling, the sequence of 4 exposures 
was duplicated for the IR observations of IC 342 and 
M33. 

We obtained the reddening E(B—V ) for each 
galaxy from the latest Galactic dust maps 
(jSchlaflv V Finkbeiner i 1201 1). and we calculated the 
extinction due to Galactic dust in each band for all of the 
galaxies in the sample using the PyR AF/STSDAS task 
Calcphot (iBushouse fe Simonl 1 994 ). We applied the 
appropriate r edde ning to an F star template spectrum 
using the iCardelli et al.l ( 1989 ) reddening law. From 
simulated T/5T/WFC3 observations, we obtained the 
magnitudes of the reddened and unreddened template 
star in each band, and determined the extinction from 
the difference between the two. All absolute magnitudes 
listed in this paper have been corrected for Galactic 
extinction and are given in the Vega magnitude system. 
The conversion between magnitudes in the Vega system 
and magnitudes in the AB system (tuab — myega) is 
1.50 mag in F275W, 1.18 mag in F336W, —0.15 mag in 
F438W, 0.42 in F814W, 0.96 in F127M, and 1.25 mag 
in F153M. 

2.3. Data Reduction 

The standard HST calibration pipeline was used for 
flat fielding, bias correction and dark current subtrac¬ 
tion. In order to reject cosmic-ray hits and bad pixels 
and to correct for the geometric distortion introduced 
by the optics of HST and WFC3, the four offset expo¬ 
sures for each filter were combined into a final image 
using the PyRAF/STSDAS task AstroDrizzle, using a 
square kernel with pixfrac = 1, exposure time weight¬ 
ing, and sky subtraction turned off. In addition, as part 
of the AstroDrizzle process, we resampled the IR im¬ 
ages onto the UVIS grid for a uniform pixel size of 0704 
in all bands. The resulting images are shown in Figure 
O Color composite images of each cluster, along with 
images of their host galaxies, are shown in Figure [3j 

3. SURFACE BRIGHTNESS PROFILE FITTING 

A 2D model of the surface brightness profile was fit 
to each i mage using GALFIT version 3.0.4 ([Peng et al.l 
120021 . I2010D . For each image, we specified the compo¬ 
nents of the model surface brightness profile to fit to the 
image (e.g., Sersic, exponential) and supplied initial es¬ 
timates for the model parameters (e.g., total magnitude, 
























Throughput 


4 


TABLE 1 

The Galaxy Sample 


Object 

D (Mpc) 

m — M Reference 

! Morphology 

u r (km s x ) 

Host Galaxy Mb 

E(B— V) 

IC 342 

3.3 ±0.3 

27.58 ±0.18 

1 

SAB(rs)cd 

31 

-23.99 

0.480 

M33 (NGC 598) 

0.9 ±0.1 

24.65 ±0.19 

2 

SA(s)cd 

-179 

-19.05 

0.036 

NGC 247 

3.4 ±0.3 

27.68 ±0.20 

2 

SAB(s)d 

156 

-18.82 

0.016 

NGC 300 

2.0 ±0.3 

26.48 ±0.28 

2 

SA(s)d 

144 

-18.04 

0.011 

NGC 2403 

3.1 ±0.2 

27.43 ±0.15 

3 

SAB(s)cd 

133 

-19.14 

0.034 

NGC 2976 

3.6 ±0.1* 

27.77 ±0.07 

4 

SAc pec 

1 

-17.74 

0.064 

NGC 3621 

7.3 ±0.2 

29.30 ±0.06 

3 

SA(s)d 

730 

-20.38 

0.068 

NGC 4244 

4.3 ±0.1* 

28.16 ±0.07 

4 

SA(s)cd 

244 

-18.96 

0.018 

NGC 4395 

4.3 ±0.4 

28.17 ±0.18 

5 

SA(s)m 

319 

-17.66 

0.015 

NGC 7793 

3.4 ±0.1 

27.68 ±0.05 

6 

SA(s)d 

227 

-18.38 

0.017 


Referen CES. — Re dshift-independent dista nces and dis t ance modu li (m — M) w e re obt ained from the fo l lowing 
Saha et ahl (120021 ): (2) IBono etahl (j2010D : (3) ISaha et all ([2006) ; (4) I Jacobs et al.l (2009); (5) IThim et al.l (2004) 


Pietrzvhski et aLI (12010 ). All distances were measured using Cepheid variables except those mark ed with an asterisk (* 


(i) 

( 6 ) 


which are TRGB distances. The E(B—V) reddening due to Galactic dust for each galaxy is from Schlaflv & Finkbeiner 

(f20TTT ). 


TABLE 2 

HST/ WFC3 Filter Properties & Exposure Times 


Filter 

Channel 

Central 

Wavelength (A) 

Band 

Width (A) 

Exp. Time 

(s) 

IC 342 Exp. 
Time (s) 

M33 Exp. 
Time (s) 

F275W 

UVIS 

2710.2 

164.51 

4 x 380 

4 x 130 

4 x 130 

F336W 

UVIS 

3354.8 

158.44 

4 x 290 

4 x 55 

4 x 55 

F438W 

UVIS 

4326.5 

197.30 

4 x 115 

4 x 35 

4 x 42 

F547M 

UVIS 

5447.4 

206.22 

4 x 95 

4 x 25 

4 x 90 

F814W 

UVIS 

8029.5 

663.33 

4 x 75 

4 x 15 

4 x 40 

F127M 

IR 

12740.0 

249.55 

4 x 60 

8x5 

8x5 

F153M 

IR 

15322.0 

378.94 

4 x 60 

8x5 

8x5 


Note. — The central wavelength listed is the “pivot wavelength”, a source-independent 
measure of the characteristic wavelength of the filter, and the band width is the “passband 
rectangular width”, the integral with respect to wavelength of the throughput across the filter 
passband divid ed by the maximum throughput, as defined in the WFC3 Instrument Handbook 
(lDresselll2012h . 



Fig. 1.— Throughput curves for the seven HST/ WFC3 filters 
taken from the PyRAF/STSDAS task Calcband. Also shown in 
gray is a model spect rum of a 1.4 Gyr old stell ar population with 
solar metallicity from Bruzual & Chariot (2003). 

position of the center of the profile). A model image 
based on these parameters is produced and then con¬ 
volved with the corresponding HST /WFC3 PSF. The 


X 2 difference between the HST image and model im¬ 
age is computed, and minimization of x 2 is achieved 
(and thus, best-fitting model parameters are found) us¬ 
ing the Levenberg-Marquardt downhill-gradient method. 
To properly weight each pixel in the y 2 calculation, we 
used the sigma image generated internally by GALFIT 
based on each input image. 

We fit 2D surface brightness profiles to the im¬ 
ages rather than fitting ID radial profiles to elliptical 
isophotes of the images, because it allowed us to more 
accurately disentangle the instrumental PSF from the 
intrinsic profile of the clusters. The clusters are quite 
compact and the fits are therefore very sensitive to the 
details of the PSF profile. HST/WFC3 PSFs have non- 
axisymmetric features that cannot be described by a ID 
radial profile. In addition, we find that some clusters 
have strongly non-axisymmetric features such as flat¬ 
tened disk components, which cannot be modeled ac¬ 
curately in ID. PSF convolution in ID is therefore not 
a suitable alternative to 2D convolution for the purpose 
of extracting the intrinsic 2D structure of the NCs from 
the HST images. 

3.1. PSF Models 

Simulated images of the HST/WFC3 PSFs in each 
filter were genera ted using version 7.4 of the TinyTim 
software package (|Kristl [19931 : iKrist et all 1201 lib which 
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Fig. 2. — Inverted grayscale images of the inner 8" X 8" of each galaxy in all seven HST/ WFC3 filters, shown with an asinh stretch. 
Images are displayed with their native orientation. 
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Fig. 3.— Color composite images of each NC and its host galaxy. For each object, the top panel shows the entire host galaxy, the middle 
panel shows the inner 12" x 12" of the galaxy , and the bottom pan el shows the inner 2" x 2" centered on the NC. Host galaxy images 
were taken from the following sources: 2MASS (Skrutskie et al. 2006, IC 342), GALEX (IGil de Paz et al.ll2007l . M33), SDSS (Baillard et al.l 
[20TT1 . NGC 2403, NGC 2976, NGC 4244, and NGC 4395), and CGS (IHo et al.N20TTl . NGC 247, NGC 300, NGC 3621, and NGC 7793). 
For the NC images, the HST/ WFC3 F275W, F547M, and F814W images were used to populate the blue, green and red color channels, 
respectively. The yellow boxes on the host galaxy images show the full 20" X 20" UVIS field of view and the scale bar denotes 1'. All of 
the images have been rotated so that north is up and are shown with an asinh stretch. 
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creates a model HST PSF based on the based on the 
instrument, detector chip, detector chip position, and 
filter used in the observations. To ensure that the model 
PSFs were processed in the same way as the HST im¬ 
ages, we produced four versions of each PSF on a subsam¬ 
pled grid with sub-pixel offsets, using the same four-point 
box dither pattern as the HST /WFC3 exposures. Each 
model PSF was convolved with the appropriate charge 
diffusion kernel in order to account for the effect of elec¬ 
trons leaking into neighboring pixels on the CCD. The 
PSFs were then combined and resampled onto a final grid 
with a pixel size of 0704 in all bands using AstroDrizzle, 
with the settings described in Section 2.3. The resulting 
model PSFs are shown in Figure HJ 

3.2. Fitting Method 

The NCs sit at the center of their host galaxies, whose 
surface brightness profiles can generally be described by 
an exponential disk model because they have either a 
weak bulge component or no bulge component at all. The 
disk of the host galaxy contributes a significant amount 
of flux in our images, particularly in the outskirts of the 
NCs, and can have a strong effect on NC fit parame¬ 
ters such as the Sersic index and effective radius. Our 
images have a small field of view and only sample the 
very inner regions of the galaxies. For example, the most 
distant galaxy in our sam ple, NGC 3621, has ap parent 
angular dimensions 12' x 7' (iGil de Paz et al.ll2007lh com¬ 
pared to the 20" x 20" UVIS field of view. Although the 
NCs sit against a non-uniform background, it is not al¬ 
ways possible to detect a gradient in the background level 
across the image due to the small field of view. In order 
to model the host galaxy surface brightness profile, we 
masked out the NC as well as any prominent dust lanes 
and attempted to fit a model consisting of an exponen¬ 
tial disk and a flat sky level to the host galaxy. If this fit 
converged, we performed the NC fit with the host galaxy 
exponential disk and sky parameters fixed at the best-fit 
values found in the previous step. We were only able to 
fit the host galaxy with an exponential disk and flat sky 
level using this method for IC 342, NGC 3621, and NGC 
4244. For the rest of the sample, the scale length and 
total magnitude were flagged as a problematic parame¬ 
ters by GALFIT, indicating that the host galaxy profile 
is too flat across the image to be modeled by an expo¬ 
nential disk profile. For M33, NGC 247, NGC 300, NGC 
2403, NGC 2976, NGC 4395 and NGC 7793, we used a 
flat sky level to model the host galaxy light and left its 
surface brightness as a free parameter when fitting the 
NC surface brightness profile. 

For each fit, we used a 401 x 401 pixel (16" x 16") 
fitting region centered on the brightest pixel in the NC. 
This fitting region size is quite large compared to the 
typical angular size of NCs in our sample, whose effec¬ 
tive radii are all well under 1". We experimented with a 
range of different fitting region sizes to quantify the effect 
of the fitting region size on the fit results and find that 
these effects are at the 10 % level for the nuclear clus¬ 
ter properties. For a single Sersic component fit with a 
flat sky level, increasing the size of the fitting region will 
lead to an increase in the brightness of the Sersic com¬ 
ponent, an increase in the effective radius, an increase in 
the Sersic index, and a decrease in the sky level. For ex¬ 
ample, when we performed a single Sersic component fit 


to the F814W image of NGC 2403 using a 150 x 150 pixel 
fitting region, we obtained a magnitude of 15.48 mag, an 
effective radius of 10.00 pixels, and a Sersic index of 1.51 
for the NC component. When we performed the same 
fit using a 300 x 300 pixel fitting region, we obtained a 
magnitude of 15.41 mag, an effective radius of 10.82 pix¬ 
els, and a Sersic index of 1.64 for the NC component. In 
order to be able to compare fit results in different filters, 
it is important to use a consistent fitting region. In order 
to optimally sample the host galaxy light, we used the 
largest fitting region size that could be applied to every 
image while still excluding bad pixels on the edges of the 
image introduced during the image combination step of 
AstroDrizzle. 

We experimented with different functions available in 
GALFIT to fit to the WFC3 images. Exponential, Gaus¬ 
sian, and de Vaucouleurs profiles are not flexible enough 
to describe the wide range of inner profil e slo p es see n 
in the NCs in our sample . Mo dified King (|Elsonl [l999b 
Nuker (|Lauer et al.l [1995) and iSersid (|1968[ ) profiles all 
provided good fits to the F814W images of M33, NGC 
300, and NGC 7793. These NCs have relatively simple 
structures and serve as instructive cases for comparison 
of different surface brightness models. In each case, the 
Sersic profile results in the smallest y 2 , although the dif¬ 
ference in x 2 between fits performed using Sersic and 
King profiles is always at the 10% level or less. The 
Nuker profile has five independent parameters and the 
King profile has four, while the Sersic profile is described 
by only three independent parameters: the total magni¬ 
tude, effective radius (half-light radius) and Sersic index. 
We used a Sersic profile to model the light profile of the 
NCs in our sample because the Sersic profile provides 
the best flexibility in describing the structure of differ¬ 
ent NCs using the smallest number of free parameters. 
To allow for elliptical shapes we left the axis ratio (5/a, 
where b and a are the lengths of the semiminor and semi¬ 
major axes, respectively) and the position angle (PA) of 
the semimajor axis as free parameters in the fits. 

With the exception of NGC 2976, NGC 4395 and NGC 
4244, a single Sersic profile was used to model each NC in 
every band. In the bands blueward of F814W, the struc¬ 
ture of the NC in NGC 2976 is difficult to fit with any 
reasonable set of analytic models due to its small-scale 
dumpiness. This NC also contains a compact blue fore¬ 
ground component to the north (see Figure [3]). For NGC 
2976, we used aperture photometry to measure the total 
magnitude of the cluster and the magnitude of the blue 
component in F275W, F336W, and F438W, and F547M. 
Since the blue component becomes faint and the struc¬ 
ture of the NC becomes smoother at longer wavelengths, 
we used GALFIT to fit a single S ersic profile to the cluster 
in F814W, F127M and F153M. SetheLal] il2008Hi found 
that the NC in NGC 4244 has two distinct morphological 
components: a red/old, compact, spheroidal structure, 
and a blue/young, extended disk stru cture, while stud¬ 
i es by lFTlippenko fc Sargent] (jl989[ ) and lFilippenko fe Hoi 
(j2003[ ) have shown that NGC 4395 hosts a NC as well as 
an unobscured AGN. Based on these studies, we modeled 
the NC in NGC 4244 using two Sersic components, and 
for NGC 4395 we used a Sersic component to model the 
NC and a point source component to model the AGN. 
For comparison to the rest of the sample, we also per¬ 
formed single Sersic component fits to the NGC 4244 and 
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Fig. 4.— Model PSFs in each filter. The full width at half maximum of each PSF is given in arcseconds. 


NGC 4395 images. A more detailed discussion of indi¬ 
vidual objects is given in Section 5. Figure [5] shows the 
results of the fits to the F814W images for each NC. 

3.3. Comparison with I SHAPE 

iGeorgiev fc Bokerl ([2014) use d the I SHAPE procedure in 
the BAOLAB software package (|Larse 3 mu to measure 
effective radii, axis ratios, and magnitudes of a sample 
of 228 NCs in late-type spiral galaxies using images from 
the JJST/WFPC2 archive. Much like GALFIT, ISHAPE 
measures structural properties by minimizing the y 2 dif¬ 
ference between an image and a PSF convolved model, 
but it performs the fit over a circular, rather than a rect¬ 
angular fitting region. In their study, the NCs were mod¬ 
eled using single component King and power-law models 
and the HST /WFPC2 model PSFs were generated using 
TinyTim. However, they did not allow for the possibility 
of an exponential disk profile to model the host galaxy 
light. The substructure seen in some of the clusters in our 
sample would not be easily detected at the distances typi¬ 
cal of NCs in their sample, which was selected using a dis¬ 
tance cut of 40 Mpc. As a test, we used ISHAPE to fit sin¬ 
gle Sersic component models to our F814W images using 
a 200 pixel fitting radius and found very good agreement 
with our GALFIT results. The mean offset between magni¬ 
tudes measured using ISHAPE and magnitudes measured 
using GALFIT (i.e., |togalfit - toishapeI) was 0.05 mag, 
while on average, the effective radii and Sersic indices 
agreed within 11% and 24%, respectively. 

3.4. ID Radial Profiles 

The “goodness of fit” is easy to visualize in a ID 
radial profile plot, although non-axisymmetric features 
of the surface brightness profiles are lost in this repre¬ 
sentation. We fit elliptical isophotes (ellipses of con¬ 
stant surface brightness) to both the HST /WFC3 im¬ 
ages and the PSF-convolved GALFIT mo dels using the 
PyRAF/STSDAS routine Ellipse (iFreudlingj [1 993]), 
which is based on methods detailed bv iJedrzeiewskil 
dm3). For each isophote, this routine computes the 
semimajor axis (SMA) of the ellipse, as well as the 
isophote surface brightness and ellipticity (e = 1 — 6/a, 
where 5/a is the axis ratio). We also performed isopho- 
tal fits to the model PSFs in order to display their radial 
surface brightness profiles. Plots of the F814W radial 
profiles of the surface brightness and ellipticity for each 
cluster are shown in Figure [6] In every case the PSF- 
convolved GALFIT models are significantly more extended 
than the PSF and the cluster effective radius exceeds the 
radius which encloses half of the flux of the F814W PSF 
(0'/063), indicating that all of the NCs in our sample are 
spatially resolved. 

3.5. Uncertainties in Cluster Parameters 


Because the data have very high S'/TV, the formal 
GALFIT measurement uncertainties in the best-fit param¬ 
eters obtained from propagating the uncertainty in elec¬ 
tron counts at each pixel are quite small. Our error bud¬ 
get is dominated by uncertainties that arise due to the 
degeneracy between the NC and the background light of 
the host galaxy light, which we model in GALFIT using ei¬ 
ther a flat sky level or the combination of a flat sky level 
and an exponential disk. Estimates of the uncertainties 
of the best-fit parameters should therefore reflect how 
sensitive these parameters are to variations in the back¬ 
ground level. In general, there is no unambiguous way to 
define a boundary between the NC and the underlying 
background light from the disk of the host galaxy. In 
addition, the background can be quite non-uniform and 
lumpy, due to the presence of resolved stars, off-nuclear 
clusters, and patchy dust lanes in the immediate vicinity 
of the NC. 

The formal uncertainty on the sky parameter returned 
by GALFIT is typically small (< 0.1 ADU), and does not 
reflect the true uncertainty resulting from the lumpy na¬ 
ture of the immediate surroundings of the NC. In order 
to more realistically estimate the uncertainty in the sky 
parameter (cr s k y ), we took the 401 x 401 pixel section of 
each image that was used in the fit, and divided it into 
eight equally sized wedges, by splitting it into quadrants 
along the vertical and horizontal axes of the image, and 
then splitting each quadrant in half at a 45° angle with 
respect to the vertical and horizontal axes. We found 
the mean sky level in each wedge while masking the NC 
and computed the standard deviation of these values, to 
obtain cr s k y . This uncertainty (typically a few ADU), 
is much larger than the formal uncertainty on the sky 
parameter returned by GALFIT. We used the mean sky 
level in this calculation, as opposed to the mode sky, 
because the mean sky is a better representation of the 
surface brightness of the host galaxy, which includes re¬ 
solved sources, such as bright stars and off-nuclear clus¬ 
ters. Using the mode sky would underestimate the host 
galaxy surface brightness by excluding the contribution 
of individual resolved stars. 

For fits where a flat sky level was used to model the 
background light from the host galaxy, the correspond¬ 
ing uncertainties in the NC best-fit parameters were es¬ 
timated by re-running GALFIT while holding the sky pa¬ 
rameter fixed at sky — a s k y (where sky is the best-fit sky 
parameter from the original fit), and again with the sky 
parameter fixed at sky + cr s k y . For fits where a combina¬ 
tion of a flat sky level and exponential disk were used to 
model the host galaxy, we re-ran the host galaxy fit as 
explained in Section 3.2, but held the sky parameter fixed 
at sky — <r s k y , and again with the sky parameter fixed 
at sky + cr s k y . We then performed the NC fit with the 
host galaxy exponential disk and sky parameters fixed 
at the best-fit values found in the previous step to find 
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Fig. 5. — Results of the F814W fits. The panels (from left to right) show the F814W image, the best-fit model convolved with the PSF, 
and the residuals (image — model). The field of view is 16 7/ x IQ" for all fits. Images are displayed with their native orientation, as in 
Figure [2] 
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Fig. 5.— Results of the fits in F814W (continued). 
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Fig. 6. — F814W radial profiles of the surface brightness and ellipticity of the GALFIT model (solid black line) and HST/ WFC3 image 
(black squares) as a function of SMA, as well as the subcomponents of the fit. These radial profiles were obtained using the elliptical 
isophote fitting routine Ellipse. The host galaxy model light profile is denoted by a blue dashed line, and the NC model light profile is 
denoted by a green long-dashed line. For reference, we include a radial profile of the F814W PSF, indicated by a red dotted line. The 
vertical gray long-dashed line shows the effective radius of the NC. The convolved NC radial profiles are significantly more extended than 
the PSF, demonstrating that the clusters are spatially resolved in each case. 


the upper and lower uncertainties in the NC best-fit pa¬ 
rameters. These uncertainties were added in quadrature 
to the formal measurement uncertainties computed by 
GALFIT in order to produce the total uncertainty on each 
fit parameter. 

The uncertainty on the distance to each galaxy is an 
additional source of systematic uncertainty that affects 
absolute magnitudes in different bands in the same way. 
In order for the uncertainties that we quote to reflect only 
our own analysis, we did not include the uncertainty on 
the distance modulus in the error budget for absolute 
magnitudes. However, interested readers can refer to the 
uncertainties on distance moduli quoted in Table [I] 

4. RESULTS 

4.1. Single Band Results 

4.1.1. Structural Properties 

The NCs in our sample have a wide range of structural 
properties. From single Sersic component fits, we find 


that F814W absolute magnitudes range from —11.20 mag 
to —15.05 mag with a mean of —12.56 mag, effective radii 
range from 1.38 to 8.28 pc with a mean of 4.27 pc, Sersic 
indices range from 1.63 to 9.73 with a mean of 3.93, and 
axis ratios range from 0.57 to 0.94 with a mean of 0.83. 
The best-fit parameters in each band are shown in Table 
02 We also include the extinction along the line of sight to 
each galaxy in this table. To help visualize the properties 
of the sample as a whole, we constructed histograms of 
the best-fit absolute magnitudes, effective radii, Sersic 
indices, and axis ratios for the single Sersic component 
fits in each band, which are shown in Figure [71 
We measure a much wider range of best-fit Sersic in¬ 
dices than in previous studies, including some very large 
values (n > 8). For the small number of NCs for 
which a measurement of the Sersic index has been per¬ 
formed, previous studies have generally fo und n ~ 1-3 
([Graham fe SpitleJ 120091: iSeth et ahll2QlQ[ ). In order to 
investigate the robustness of these high concentrations, 
we used GALFIT to produce noise-free Sersic models with 
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Fig. 6.— F814W radial profile plots (continued). NGC 4244 is fit with 2 Sersic components, which are indicated by a green long-dashed 
line and an orange dash-dotted line. 


an effective radius of 0'/23 (which is the mean value for 
our sample) and Sersic indices ranging from 0.5 to 10, all 
which were convolved with the F547M PSF. We show 
their radial surface brightness profiles, along with ra¬ 
dial profiles of the corresponding intrinsic Sersic models 
(which have not been convolved with the PSF) in Figure 
[8l The large values of the Sersic index that we mea¬ 
sure are determined more by the behavior of the profile 
in the wings than in the core, because the PSF blurs 
out information about the slope of the inner profile. For 
radii inside the half width at half maximum of the PSF, 
the profiles for n > 6 are not distinguishable. However, 
these profiles are strongly divergent at radii outside of 
the core of the PSF. Therefore, the large Sersic indices 
that are seen in some objects serve as useful indicators 
of the shape of the NC profile outside of the core of the 
PSF. We conclude that NCs have a wide range of concen¬ 
trations as indicated by the wide range of Sersic indices, 
and that Sersic indices above 6 do not indicate differ¬ 
ences in core concentrations, but instead are dependent 
on the wings of the profile. 

4.1.2. Size vs. Luminosity 


We show the F547M effective radius plotted against 
the F547M absolute magnitude for the NCs in our sam¬ 
ple (Figure [9|) for comparison with the NCs and the 
derive d size-luminosity relation from iGeorgiev fe Bokerl 
(l20li ). The magnitudes that they measured in 
the WFPC2 F606W and F555W filters were con¬ 
verted to T/5T/WFC3 F547M magnitudes using the 
PyRAF/STSDAS task Calcphot, assuming a 5 Gyr, so¬ 
lar metallicity simple stellar population (SSP) for the 
spectral shape (this is the same spectrum that they used 
to convert between HST /WFPC2 magnitudes and John¬ 
son magnitudes). For NGC 2976, we show the F547M ab¬ 
solute magnitude and the F81 4W effective radius . Most 
of our clusters fall below the IGeorgiev fc Bokerl ( 2014 ) 
best-fit size-luminosity relation, which may be in part 
due to our small sample size. In addition, because they 
used a large distance cut (40 Mpc), 11% of the NCs in 
their sample were not spatially resolved and thus were 
not used to derive their size-luminosity relation. As a 
result, the subset of spatially resolved NCs in their sam¬ 
ple may be somewhat biased towards larger objects. It 
is also worth noting that the NC in IC 342 is exception¬ 
ally bright and compact. This NC is by far the bright- 
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Fig. 6 .— F814W radial profile plots (continued). NGC 4395 is fit with a Sersic component (indicated by a green long-dashed line) 
and a point source component (orange dash-dotted line), and the vertical gray long-dashed line denotes the effective radius of the Sersic 
component of this two component fit. We do not include an additional PSF for reference in this plot because the point source is included 
a fit component. 


TABLE 3 

NC BEST-FIT PARAMETERS 



F275W 

F336W 

F438W 

F547M 

F814W 

F127M 

F153M 

IC 342 A (mag) 

2.94 

2.44 

2.00 

1.51 

0.90 

0.41 

0.30 

m (mag) 

1 7 9^+0.25 

l/.ZO_Q Qg 

16.21+23 

15.99+°® 
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11 14+ 0 - 06 
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M (mag) 
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lo.ZjI _q 08 

— 13.80+i° 
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R e (pc) 
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i -^ Z - 0.12 

1 qo+ 0-48 
1.58—o i 6 

1 41+0.38 

1 42+ 0 - 33 

1 qo+ 0.19 
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4 -^ U -0.09 

1 40+ 0 - 20 
i>4:U -0.09 

n 

i-03lS:g 

l-94+gj 

2-56+®® 

2 Q7 + 0 - 78 
'-0.36 

4 7Q + 0 - 48 

1 y — 0.40 

2 47 + 0 - 73 
z -^'-0.32 

2 4Q + 0 - 69 
Z '^ y -0.32 

b/a 

0.98+“ 

0.84+“ 

0.89+“ 

0 Q 2+ 0 - 02 
u - yz - 0.01 

n 04 + 0.01 
u - y 4 : - 0.01 

O 91 +o.oi 
u - yi - 0.01 

0 Q2 + 0 - 01 
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Note. — Best-fit parameters and their corresponding uncertainties for each NC. We also include the extinction in each 
band ( A ) along the line of sight to each galaxy. Some axis ratios have very small uncertainties, which round to 0.00. We list 
an uncertainty of 0.01 in these cases. Best-fit parameters for all 10 NCs are available in the online version of the table. 


est in our sample, and has a F547M absol ute magnit ude 
bright er than 95% of the NCs in the iGeorgiev fc Bokerl 
( 2014 ) sample. 

4.1.3. NC Magnitude vs. Host Galaxy Magnitude 

IBoker et al.l ( 2004 ) found that the luminosities of nu¬ 
clear clusters in late-type spiral galaxies generally cor¬ 
relate with the total luminosities of their host galax¬ 
ies. In particular, the /-band luminosities of NCs show 
a strong correlation with the 5-band luminosity of the 
host, and a weak correlation with the far-IR luminosity 
of the host, which suggests that NC masses are more 
governed by the mass of the host galaxy rather than 
the total star formation rate of the host galaxy. In Fig¬ 
ure m we plot the /-band absolute magnitudes of the 
NCs in our sample against the 5-band absolute magni¬ 
tudes of their host galaxies, wh ich were taken from NED . 
We also include the NCs from Georgiev & Boker (2014) 
and show the best-fit relation from IBoker et al.l ([2004D . 
As described above, HST magnitudes were converted to 
Johnson magnitudes using Calcphot, assuming a 5 Gyr, 
solar metallicity SSP for the spectral shape, and extinc¬ 
tions in the Johnson 5 and / bands were determined 


using the method described in Section 2.2. Most of the 
objects in our sample lie close to the best-fit relation 
and have luminosities that are quite typical for NCs in 
late-type spiral galaxies. 

4.2. Panchromatic Results 
4.2.1. Color Gradients 

We find a strong wavelength dependence of the the ef¬ 
fective radius in the majority of the clusters, indicating 
radial color gradients. In Figure [HI we plot the best fit 
effective radius of each NC against wavelength, showing 
the results of the single Sersic component fits for NGC 
4244 and NGC 4395. We performed a formal power-law 
fit to each curve shown in Figure [TT] to compute their 
slopes [dlog(5 e /pc)/dlog(A/A)] and corresponding un¬ 
certainties. For IC 342, NGC 2976, NGC 2403, and NGC 
4244, the best-fit slope is consistent with zero at the 1 
a confidence level. Positive slopes were found for M33, 
NGC 247, NGC 300, and NGC 4395 (best-fit slopes of 
0.24 ± 0.06, 1.00 ± 0.13, 0.53 d= 0.02, and 1.00 ± 0.01, re¬ 
spectively). This suggests the presence of a young pop¬ 
ulation that is more concentrated than the bulk of the 
stars in the NC, or in the case of NGC 4395, indicates the 
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Fig. 7.— Histograms of the best-fit parameters in each band. For NGC 4244 and NGC 4395, the results of the single Sersic component 
fits are shown. From left to right, each row contains (for a given band) histograms of absolute magnitudes, effective radii, Sersic indices, 
and axis ratios. The median value for each parameter is denoted by a vertical dashed line. 
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Fig. 8 .— Sersic models for R e = O'/23, with and without convo¬ 
lution with the F547M PSF. Solid lines denote Sersic models that 
have been convolved with the F547M PSF, while the dotted lines 
denote intrinsic Sersic models which have not been convolved with 
the PSF. Convolved models are normalized to the same surface 
brightness at R = 0, while unconvolved models are normalized at 
R = 0723. The black vertical long-dashed line indicates the half 
width at half maximum of the F547M PSF. The black horizon¬ 
tal dash-dotted line indicates the ratio between the central surface 
brightness of the host galaxy and the central surface brightness of 
the NC in F814W, averaged over all ten objects in the sample. 



Fig. 9.— The effective radius measured in the F547M filter vs. 
the F547M absolute magnitude. Black triangles denote NCs from 
this work. The black “x” indicates the single Sersic component fit 
to NGC 4395 (that is, the effective radius and magnitude of the 
NC plus AGN). T h e sma ll light blue diamonds indicate NCs from 
IGeorgiev fc BdkeH (|2014T ) and the gray dashed line indicates the 
best-fit V band size-luminosity relation derived from their data 
set. 


Fig. 10.— /-band absolute magnitudes of NCs vs. B-band abso¬ 
lute magnitudes of their host galaxies. Black triangles denote NCs 
from this work, and sma ll light blue diamonds indicate NCs from 
IGeorgiev &; Bdk erjj2014lj. The ? ^ray dashed line indicates the best 
fit relation from lBoker et al. (2004). 

presence of the AGN. These results are consistent with 
results from the recent H ST /WFPC2 archival study by 
IGeorgiev fe Bokerl (j2014D . which found that on average, 
NCs have larger effective radii when measured in red¬ 
der filters. However, we find a trend of decreasing ef¬ 
fective radius with wavelength in NGC 3621 and NGC 
7793 (best-fit slopes of 0.95 =b 0.07 and —0.30 ± 0.10, re¬ 
spectively) , which may indicate recent circumnuclear star 
formation. It should be noted that the larger extent of 
NGC 3621 in UV bands may be due to UV photons from 
the central engine of the AGN scattering off narrow-line 
region clouds rather than a radial gradient in the stellar 
populations. These color gradients will be explored in 
more detail in future work, where we will perform spa¬ 
tially resolved stellar population modeling of the NCs in 
this sample. 

4.2.2. Increasing Roundness With Wavelength 

In a majority of the clusters, we also find a trend of 
increasing roundness with wavelength. In Figure [12j 
we plot best-fit NC axis ratios vs. wavelength, show¬ 
ing the results of the single Sersic component fits for 
NGC 4244 and NGC 4395. We also include the axis ra¬ 
tios of the PSF in each band, all of which are nearly 
unity (> 0.98). We performed a formal power-law fit to 
each curve shown in Figure [12] to obtain best-fit slopes 
[dlog(b/a)/dlog(A/A)] and their corresponding uncer¬ 
tainties. For every NC in the sample, this slope is not 
consistent with zero at the 1 a confidence level. IC 342 
and NGC 7793 have slightly negative slopes, while every 
other NC has a positive slope, indicating a monotoni- 
cally increasing axis ratio with wavelength. This suggests 
that the young stars in NCs often form in flattened disks, 
regardless of whether these disks are more extended or 
compact than the older stellar populations in the cluster. 
If a cluster contains a young population of stars in a flat¬ 
tened disk structure in addition to an older population in 
a spheroidal structure, the 2D axis ratio should be closer 
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Fig. 11.— Effective radius of the NCs in each filter. For NGC 
4244 and NGC 4395, the results of the single Sersic component fits 
are shown. 



Filter Central Wavelength (A) 

Fig. 12. — Axis ratios of the NCs in each filter. For NGC 4244 
and NGC 4395, the results of the single Sersic component fits are 
shown. The gray x’s denote the axis ratio of the PSF in each filter. 

to unity when measured in redder bands, provided that 
the disk component is not face-on. On the other hand, if 
the young and old populations have identical morpholo¬ 
gies, the axis ratio should not exhibit any trends with 
wavelength, but would not necessarily be equal to unity. 
Thus, if the young stars in NCs typically form in disks, 
while the older population is in a spheroidal distribution, 
we would expect objects containing a young stellar disk 
to exhibit a monotonic increase in axis ratio with wave¬ 
length. In NGC 4244, for example, the cluster is clearly 
being viewed edge-on, and the best-fit axis ratio from 
the single Sersic component fits is 0.29to.oi i n F275W 
and 0.44lo oi i n F15 3M. T his is c onsistent with a pre¬ 
vious study by iSeth et ah l ( 2008b ), which showed that 
the NC contains morphologically distinct old and young 
populations. 

If the half-light radius of the cluster is significantly 


smaller than that of the PSF, the intrinsic axis ratio of 
the cluster will not be recovered in the fit. In the limit 
of spatially unresolved sources, the best-fit axis ratio will 
be completely determined by the axis ratio of the PSF, 
which is always near unity. For our sample, the clusters 
are sufficiently well resolved that we were able to recover 
information about their axis ratios. With the exception 
of IC 342, all the clusters have effective radii larger than 
the PSF half-light radius in all bands. This suggests the 
variations in size and flattening with wavelength that we 
are observing are robust, indicating a difference in the 
spatial distribution of young and old stellar populations. 

The flattening preferentially seen in the bluest bands 
indicates that the young stars form in disks, possibly 
fueled by the accretion of low angular momentum gas 
ont o the center of the g alaxy. However, simulations 
by lHartmann et al.l ( 20111 ) demonstrated that the star 
cluster migratory/merger scenario can also produce NCs 
with young, blue disks, such as in NGC 4244. These 
simulations also showed that the star cluster migra¬ 
tory/merger and gas infall formation scenarios predict 
different kinematics in the resulting NC. Spatially re¬ 
solved measurements of the kinematics of the NCs are 
needed to distinguish between these formation scenarios. 

In Figure [T3l we plot dlog(6/a)/dlog(A/A) (top panel) 
and F814W axis ratios (bottom panel) of the NCs against 
the axis ratios of the host galaxy disk, w hich were calcu¬ 
lated using 2MASS A-band images (jJarrett et al.ll2QQ3 ). 
The galaxies in our sample with the smallest axis ratios 
contain the NCs that have the strongest trend of increas¬ 
ing roundness with wavelength and are the most flattened 
(this can also be seen in Figure [3]). The HST/WFC3 im¬ 
ages of the NCs in NGC 2976, NGC 3621, and NGC 4244 
clearly show the presence of a disk which has the same 
orientation as the disk of the host galaxy. In each of these 
cases, the PA of the NC as measured by GALFIT is nearly 
aligned with the PA of the host galaxy (APA < 20°). 
Conversely, none of the galaxies in our sample which are 
close to face-on contain flattened NCs. We conclude that 
blue disks which lie in the plane of the host galaxy are 
a common feature of NCs in late-type galaxies, but are 
difficult to detect in galaxies that are close to face-on. 

4.3. Stellar Populations 

4.3.1. Nuclear Cluster SEDs 

For each cluster, we measured a total apparent magni¬ 
tude using either GALFIT or aperture photometry in all 
seven bands, allowing us to construct a spectral energy 
distribution (SED) spanning the near-UV to the near- 
IR in wavelength. We corrected for Galactic extin ction 
using dust maps from iSchlaflv fe Finkbeinerl ( 2011 ), as 
explained in Section 2. Magnitudes were converted to 
flux units (erg s -1 cm -2 A -1 ) using Calcphot, then fi¬ 
nally converted to luminosities using the distances given 
in Table [TJ The resulting SEDs are shown in Figure HU 
Color-color diagrams (see below) provide a way to vi¬ 
sualize the comparison between the observed SEDs and 
model spectra. Fits of multi-component stellar popula¬ 
tion models to the NC SEDs are needed to more precisely 
characterize the stellar populations present in each NC, 
and will be described in a future paper. 

4.3.2. Nuclear Cluster Colors 
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Fig. 13.— Top panel: cHog(6/a)/<i log(A/A) of each NC vs. the 
axis ratio of the host galaxy. These slopes [d\og(b/a)/d log (A/A)] 
were calculated by performing a formal power-law fit to the curves 
shown in Figure fl2] Bottom panel: axis ratios of the NCs measured 
in F814W vs. the axis ratios of their host galaxies. 


In order to characterize their stellar populations, we 
computed NC colors from our extinction-corrected abso¬ 
lute magnitudes and compared them t o co lor s of sing le¬ 
age si mple stellar populations from iBruzual &; Chariotl 
( 2003 ) for metallicity Z = 0.008, 0.02, and 0.05. We 
created two color-color diagrams, one spanning UV to 
the IR, and one restricted to the optical bands, which 
are given in Figures [15] and [16] respectively. For NGC 
4395, we show colors of the NC (Sersic component, de¬ 
noted by a black triangle) as well as the combined colors 
of the NC and AGN (Sersic plus point source compo¬ 
nent, denoted by a black “x”). We also i nclude colors 
of the NCs from iGeorgiev fc Bokerl (j2014l ) on the op¬ 
tical only color-color plot. The magnitudes that they 
measured in the WFPC2 F450W, F606W, and F814W 
filters were converted to magnitudes in the 7/5T/WFC3 



MA) 


Fig. 14. — SEDs of the 10 NCs in our sample. For NGC 4244 and 
NGC 4395, we show the composite SED of the two components of 
the fit. 

F438W, F547M, and F814W filters using Calcphot, as¬ 
suming a 5 Gyr, solar metallicity SSP for the spectral 
shape. In this plot, our NCs fall close to the SSP tracks 
and are broadly consistent w ith the colors of the NCs 
from IGeorgiev fc Bokerl ([2Q14J) . 

Most of the NCs in our sample he at locations dis¬ 
tinctly offset from the SSP tracks in F336W — F438W 
vs. F547M — F127M color-color space, confirming that a 
mix of stellar populations with different ages is present in 
many of the NCs. This is expected, given that previous 
studies have shown that the star formation histories of 
NCs a re generally not well described by a singl e starburst 
event (|Walcher et al.l 120061: iRossa et al.l 120061) . In order 
to inve stigate the mix of ages, w e took linear combina¬ 
tions of lBruzual fe Chariotl (2 0031 ) SSP spectra with ages 
11 Gyr and 25 Myr and computed their colors in order to 
construct a “mixing curve”, which traces the color of the 
total population as the fractional contribution from the 
11 Gyr population to the total luminosity is adjusted 
from 1.0 to 0.0 and the contribution from the 25 Myr 
population to the total luminosity is adjusted from 0.0 
to 1.0, in steps of 0.1. A mixing curve for a combination 
of a 11 Gyr and 290 Myr SSP is also included in the 
UV-optical-IR color-color plot. The typical age mix for 
the NCs in our sample is a combination of an old (>1.4 
Gyr) and a young (100 - 300 Myr) population. The mix 
of ages is in the NCs is obvious when considering colors 
that span the UV to the IR, but is not apparent from 
the optical colors only. These color-co lor plots reinforce 
the conclusions of lde Meulenaer et al 1 (l20ll . who high¬ 
lighted the importance of broad wavelength coverage in 
stellar population studies, and demonstrated how degen¬ 
eracies between age, metallicity and extinction, which 
limit stellar population studies that utilize optical bands 
only, can be broken with full UV to IR coverage. 

5. COMMENTS ON INDIVIDUAL OBJECTS 
5.1. IC 342 

IC 342 is a face-on giant Scd spiral galax y locat ed in 
the Maffei Group at a distance of 3.3 Mpc (jSaha et al.l 
M). In a dynamical study of IC 342 using CO band- 
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Fig. 15.— F336W - F438W vs. F547M - F127M colors. Black 
triangles denote NCs from this work, and the off-nuclear clusters 
in IC 342 are shown as black squares. For NGC 4395, we show 
colors of the NC as well as the combined colors of the NC and 
AGN (denote d by a black “x”). We in clude evolutionary tracks 
of SSPs from fBruzual & "C hariot (2003) with Z — 0.008 (denoted 
by the red dashed line), Z = 0.02 (purple short-dashed line), and 
Z — 0 .05 (orange dash - dot line). Mixing curves for a combination 
of IBruzual fe Charlotl ( 2003 ) SSPs with age 25 Myr and 11 Gyr 
are denoted by the solid gray green lines, and mixing curves for a 
combination of 290 Myr and 11 Gyr SSPs are denoted by the long- 
dashed gray green lines. We include mixing curves for metallicity 
Z — 0.008, 0.02, and 0.05. The arrow indicates a reddening vector 
of Ay = 0.2 mag. The error bars in the lower right corner show 
the mean uncertainty in the colors on each axis (averaged over all 
NCs in the sample). 



F547M-F814W (mag) 

Fig. 16.— F438W - F814W vs. F547M - F814W colors. We 
include NC colors from IGeorgiev & Bokeii (2014). denoted by small 
gold diamonds. We did not include mixing curves in this plot 
because the model tracks are nearly straight lines in this color- 
color space, especially at older ages. 


head spectra and HST/ WFPC2 images, iBdker et, al.1 
measured an integrated stellar velocity dispersion 
of 33 km s -1 for the NC, estimated the total stellar mass 
to be 6 x 10 6 M 0 , and placed an upper limit on the mass 
of any central BH of 5 x 10 5 M 0 . Due to its proximity to 
the Galactic plane, there is significant foreground extinc¬ 
tion, particularly in the UV bands (Ap275W= 2.94 mag). 
There are dust lanes to the east and west of the main 
cluster that we masked out during the fits. We were able 
to fit a model consisting of an exponential disk and a 
flat sky level to the host galaxy in every band, and the 
nuclear cluster was fit usin g a sin gle Sersic component in 
every band. iBoker et al.l ( 1999 ) measured an apparent 
U-band magnitude of 15.19 mag for the cluster, which 
corresponds to an extinction-corrected F547M absolute 
magnitude of —13.91 mag, which is somewhat discrepant 
with our measurement of —14.07^0.mag. However, the 
HST /WFC2 images that they used were saturated in the 
core of the NC, which may explain the discrepancy. The 
NC in IC 342 is one of the most compact and by far the 
most luminous in our sample, with an /-band absolute 
magnitude Mp 8 i 4 w= — 15.051 q.o 8 ma g- The effective ra¬ 
dius is nearly constant with wavelength, with a best-fit 
value of 1.38lo i 3 pc in F814W. 

There are several off-nuclear clusters in the central re¬ 
gions of IC 342. The brightest two are located 2'/11 (pro¬ 
jected distance 33.55 pc) to the north and F/92 (30.53 
pc) to the south of the NC, respectively. Although other 
galaxies in our sample contain off-nuclear clusters, these 
two were included in the color-color diagrams in Section 
4 because they are by far the most luminous in our sam¬ 
ple, and actually outshine the main NC in UV bands. 
We measured their magnitudes using aperture photom¬ 
etry over a circular aperture (and applied an aperture 
correction to account for flux lost outside the aperture 
due to PSF blurring), because GALFIT could not converge 
on a good fit without the use of constraints. These off- 
nuclear clusters have extremely blue colors, consistent 
with 5 Myr SSPs, as shown on the color-color diagrams 
shown in Figures [T5l and flGl We will carry out a more de¬ 
tailed examination, including an analysis of stellar pop¬ 
ulation ages and stellar masses, of the complete set of 
off-nuclear clusters in our sample in future work. 

5.2. MSS 

Located at a distance of 0.9 Mpc (|Bono etld1l2010f h 
M33 is the only member o f the Local Group in our sam¬ 
ple. iGebhardt et al.l (|200lD measured an integrated ve¬ 
locity dispersion of 20 km s -1 for the cluster, and placed 
a very tight upper limit on the mass of an y central BH 
of 150 0 M 0 . Using the same HST data, Merritt et al.l 
(2001) placed an upper limit of 30 0 0 M 0 on the mass of 
the central BH. iKormendv et al ] <l20lQh estimated the 
total stellar mass in the cluster to be 10 6 M 0 . The 
nucleus of M 33 coincides with th e ultraluminous X-ray 
source X -8 (|Trinchieri et al.l [l988). which is believed to 
be a stellar-mas s BH accreting at a super-Eddington rate 
(|Foschini et al.ll2004h . 

In all bands, we fit a flat sky level to the host galaxy 
and a single Sersic component to the nuclear cluster. 
The NC in M33 is one of the most compact in the 
sample. We find a steady increase in the effective ra¬ 
dius with filter central wavelength from 0.69 ^q 'to P c in 
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F275W to 1.78^o'.i3 P c i n F153M, which indicates that 
the center of the cluster is bluer than the wings. This 
is consisten t with a /7ST/WFPC2 imaging study by 
lLauer et~ahl (j1998t h which found a decrease in the V—I 
col o r toward s th e cent er of the cluster, and a study by 
iKormendv fc McClure (1993 ) which found a similar gra¬ 
dient in the B—R color. 

5.3. NGC 247 

NGC 247 is a member of the Sculptor Group, at a dis¬ 
tance 3.4 Mpc dBono et al .l 201 0 ). We fit a flat sky level 
to the host galaxy and a single Sersic component to the 
NC in all bands. We find evidence for a color gradient 
within this cluster as well. The effective radius is small¬ 
est in the UV and increases towards the IR, with a best- 
fit value of 0 . 791 q’o 7 i n F275W and 3.491 F153M. 

The best-fit axis ratio is 0.55lo‘oi i n F275W and mono- 

tonically increases to 0.891q‘.oi i n F153M, which suggests 
that the young stars formed in a disk. 

5.4. NGC 300 

At a distance of 2.0 Mpc ( Bono et al . 201 0 ). NGC 300 
lies in the field between theLocal Group and the Sculptor 
Group. IWalcher et al.l ( 2005 ) measured a stellar velocity 
dispersion of 13.3 km s -1 and a mass of 10 6 02 M 0 for the 
NC. We fit a flat sky level to the host galaxy and a sin¬ 
gle Sersic component to th e NC in all bands. NGC 300 
was included in the iBoker etldl (2 0021 ) //ST/WFPC 2 
NC imaging survey, where an /-band absolute magni¬ 
tude of —11.43 zb 0.40 mag was measured for the NC. 
Correcting for the slightly different distances used, this 
corresponds to a WFC3 F814W absolute magnitude of 
— 11.08±0.40 mag, which is consistent with our measure¬ 
ment of — 11 . 20 lo’o 2 ma g- The best-fit effective radius 
increases monotonically from 1.52 +q’ 03 pc in F275W to 

3.12lo i 6 P c i n F153M, indicating a color gradient within 
the cluster. The Sersic index is nearly constant across all 
filter bands, with a value of 2.271 q'o8 F814W. 

5.5. NGC 2403 

NGC 2403 resides in th e outsk irt s of the M81 Group at 
a distance of 3.1 Mpc ([Saha et al.ll2006[) . In addition to 
the NC, NGC 2403 hosts a highly variable central X-ray 
source, likely to be an X-ray bin ary sy stem with a stellar 
mass BH accretor ([Yukita et aLll2QQ~7h . We fit a flat sky 
level to the host galaxy and a single Sersic component to 
the NC in every band. Dust lanes southwest of the NC 
were masked out during the fits. With an effective radius 
of 6.13 pc and an absolute magnitude of — 11.62^q in 
F814W, the cluster has a low surface brightness, lead¬ 
ing to relatively large uncertainties in the fit parameters. 
The effective radius of the cluster is generally larger when 
observed in redder filters, with a value of in 3.571q 90 P c 
F275W and 1.2^^ pc in F153M. 

5.6. NGC 2976 

NGC 2 976 is located at a distance of 3.6 Mpc 
( Jacobs et al.l 2009 ) in the M81 group. In F814W, 
F153M, and F153M, we were able to fit a single Sersic 
component to the NC while fitting a flat sky level to the 


host galaxy. In these bands, the cluster is highly flat¬ 
tened and the PA is nearly aligned with that of the host 
galaxy disk ( b/a = 0.57 and APA = 13° in F814W). 
In the bands blueward of F814W, the structure of the 
NC is difficult to fit with any reasonable set of analytic 
models due to its small-scale dumpiness. In these bands, 
the cluster can be resolved into two distinct structures: 
a compact blue component to the north and a more ex¬ 
tended, irregularly shaped red component to the south 
(see Figure [3]). The contrast between these components 
diminishes at larger wavelengths, and they completely 
blend together in the F814W, F127M, and F153M im¬ 
ages (as can be seen in Figure ET In order to obtain 
the total cluster magnitude in F275W, F336W, F438W, 
and F547M, we performed aperture photometry using an 
elliptical aperture centered on the centroid of the clus¬ 
ter. We used the axis ratio and position angle of the 
best-fit Sersic model in F814W to define the position 
angle and axis ratio of the elliptical aperture. The semi¬ 
major axis of the aperture we used was 9 x R e ,F8i4W 
(where R e ,F8i4W is the best-fit effective radius for the 
cluster in F814W), while the inner and outer radii of the 
sky annulus were 9 x R e ,F 8 i 4 W and 10 x R e ,F814W, re¬ 
spectively. These aperture and sky annulus sizes were 
found to give the best agreement with the total cluster 
magnitude obtained from GALFIT in F814W, F127M, and 
F153M. We calculated a mean difference of 0.10 mag be¬ 
tween the total cluster magnitude obtained from GALFIT 
and from aperture photometry in F814W, F127M and 
F153M, and used this to define the measurement uncer¬ 
tainty. From aperture photometry, we measured the fol¬ 
lowing extinction-corrected absolute magnitudes for the 
cluster: Mp275W = —9.10 mag, Mp336W = —9.43 mag, 
A/f 438 W = —9.13 mag, A/f 547 M = —9.97 mag. 

We also measured the magnitude of the blue compo¬ 
nent in F275W, F336W, F438W, and F547M by per¬ 
forming aperture photometry over a circular aperture. In 
each band, we performed the measurement 12 times using 
aperture radii in the range 2-5 pixels, inner sky radii in 
the range 5-20 pixels and outer sky radii in the range 10- 
30 pixels, applying the appropriate aperture correction 
each time. We calculated the median between the largest 
and smallest magnitude obtained from the 12 trials, and 
used this range to define the associated measurement 
uncertainty. We obtained the following absolute mag¬ 
nitudes for the blue component: Mp 275 W = —8.49 ±0.09 
mag, A/f 336 W = —8.38 ± 0.22 mag, A/f 438 W —7.35 ± 

0.50 mag, M^m 7 M = —7.98 ± 0.78 mag. For compari¬ 
son, the brightest Wolf-Rayet and O stars in the R136 
cluster in the LMC h ave V -band absolute m agnitudes as 
bright as -7.33 mag (jCrowther et al.l [2010). Because of 
the large uncertainties in our measurements, we cannot 
rule out the possibility of the blue component being an 
extremely luminous, hot star, although it may also be 
a very young, compact star cluster. We are unable to 
determine whether this object is associated with the NC 
or being viewed in projection. 

5.7. NGC 3621 

NGC 3621 is the most distant galaxy in our sample, 
at 7.3 Mpc ([Saha et al.l 1 20061 ). The NC has a velocity 
dispersion of 43 km s -1 , one of the largest ever measured 
for any NC in a late-type spiral galaxy, which sets a 
fairly high upper limit on the mass of any central BH of 
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3 x 10 6 M 0 flBarth et al.ll2009[) . iSatvapal et al.l ([20071 ) 
discovered a faint AGN in the NC from mid-IR [Ne V] 
line emission and based on estimates of its bolometric 
luminosity, set a lower limit o n t h e mass of the central 
BH of 4 x 10 3 M 0 . iGliozzi et al.l ( 2009 ) found further 
evidence for an AGN in NGC 3621 with the detection of 
a weak X-ray point source coincident with the nucleus. 

We successfully fit an exponential profile plus a flat 
sky level to the host galaxy in every band, and we fit a 
single Sersic component to the NC in all seven bands. 
Although NGC 3621 contains an AGN, adding a point 
source component did not improve the quality of the 
fits, consistent with expectations for a highly obscured 
AGN. Dust lanes northeast, southeast, and nort hwest of 
the NC were masked out during the fits. iBarth et al.l 
( 2009 ) measured a J-band apparent magnitude of 14.14 
mag for the NC from archival NICMOS images, which 
corresponds to a F127M absolute magnitude of —15.24 
mag, in very good agreement with our measurement of 
— 15.25to.i3 mag. The NC has a very large Sersic in¬ 
dex (n > 7) in all bands, with a best-fit value as high 
as 9.73^‘67 in F127M. In the UV, the cluster is highly 
elongated along the north-south axis, with an axis ratio 
( b/a ) of 0.48to.oi• This elongated UV emission can also 
be seen in the color composite image shown in Figure [3j 
The axis ratio monotonically increases with wavelength 
to 0.93lo‘oi F153M. The cluster is also more extended 

in the UV, with an effective radius of 17 . 27^333 pc in 

F275W and in 6.15^0.72 F127M. The greater north-south 
elongation and larger extent of the cluster in the UV 
could potentially be explained by UV emission from the 
central engine being obscured by an edge-on torus ori¬ 
ented along the east-west axis and scattering off narrow- 
line region clouds to the north and the south. If so, the 
extended UV emission would be linearly polarized. Al¬ 
ternatively, the extended UV emission could be explained 
by extended star formation in the outskirts of the clus¬ 
ter. In F275W, the PA of the semimajor axis of the 
NC is nearly aligned (APA = 19°) with the PA of the 
host galaxy m easured from the 2 MASS A-band image 
(jJarrett et al.l 2003), which is expected if the NC primar¬ 
ily grows through gas accretion rather than star cluster 
mergers (jSeth et all 120061: [Hartmann et al.ll20llt ). How¬ 
ever, there is no correlation b etween t he ori enta tion of 
AGNs and their host galaxies ([Fischer et al.ll2013f ). so we 
cannot determine the nature of the extended UV emis¬ 
sion from our UV images alone. UV spectroscopy and/or 
polarimetry can help to distinguish between young stars 
and scattered AGN emission. 

5.8. NGC 4244 

Located at a distance of 4.3 Mpc ([Jacobs et al.l [2009 ) 
in the M94 Group, NGC 4244 is the only galaxy in our 
sample with an edge-on orientation. In a recent Gemini 
Near-Infrared Integral Field Spectrograph (NIFS) study 
of this NC, iSeth et all ( 2008ld ) discovered rotation of 30 
km s -1 in the central 10 pc of the cluster and found that 
the NC contains two distinct populations: an old, com¬ 
pact spheroidal component and a young, extended disk 
component. In each band, we fit an edge-on exponential 
disk plus a flat sky level to the host galaxy, while the NC 
was modeled with two Sersic components. In F814W, the 


flattened component has an effective radius of 6.12lQ 3 g 
pc, a Sersic index of 0.521^045 an d an ax ^ s ra ti° of 
0.20to 02 > while the rounder component has an effective 
radius of 3.53to.3i pc, a Sersic index of 1.7lto.20, an d an 
axis ratio of 0.67to‘.oi- SEDs of the two components, as 
well as plots of their effective radii, Sersic indices, and 
axis ratios as a function of wavelength are shown in Fig¬ 
ure EH We also fit a single Sersic component to the NC 
in each band for comparison with other clusters. The sin¬ 
gle Sersic component effective radius is nearly constant 
with wavelength, with a best-fit value of 4.771 q!o2 P c i n 
F275W and 4.98^o’.oi F153M. From the single Sersic 

component fits, we find that the axis ratio monotonically 
increases from 0.291 q!oi F275W to 0.44tpm in F153M, 

and the semimajor axis is nearly aligned with that of the 
host galaxy in every band (APA = 12° - 13°). The 
alignment of the NC and host galaxy semimajor axes is 
consistent with previous studies of NGC 424 4 Seth et al.l 

120061 : lHartmann et al.ll2011l : iDe Lorenzi et al.ll2013[ ). and 

suggests that gas accretion has played a significant role 
in the formation of the NC. 

5.9. NGC 4395 

NGC 4395 is ori ente d face-on at a distance of 4.3 
Mpc (|Thim et al.l l2004fl in the M94 Group. The ex¬ 
istence of a NC in NGC 4395 w as explicitly demon¬ 
strated bv iFilippenko fe Hoi (|2QQ3h . from GALFIT anal¬ 
ysis of a HST/ WFPC2 /-band image. The discovery of 
faint broad hydrogen and helium lines in its optical spec¬ 
trum revealed that NGC 4395 co ntains one of the least 
luminous Type 1 AGNs known (|Filippenko fe Sargentl 
119891 ). A HST /STIS reverberation-mapping program by 
iPeterson et al.l (|2QQ5f ) yielded an estimate of 3.6 x 10 5 
Mq for the mass of the central BH, which implies that 
the AGN is very under luminous and is accreting at just 
1.2 x 10 -3 times the Eddington rate. In a r ece nt pho- 
tometric reverberation mapping campaign, lEdri et ahl 
( 2012 4 measured an Ha lag of just 3.6 hours, correspond¬ 
ing to a BH mass of 5 x 10 4 M 0 (the large discrepancy in 
the BH mass estimates from this study and Peterson et 
al. is in part due to the different assumed values of the 
scaling factor in the BH mass calculation, which depends 
on the unkno wn geometry an d kin ematics of the bread¬ 
line region). [Minezaki et al.l (2006) detected J and H 
band flux variations over a single night of monitoring of 
the AGN, in addition to large flux variations in optical 
to near-IR bands on timescales of days to months. 

We used a flat sky level to model the host galaxy, a sin¬ 
gle Sersic component for the NC, and included a point 
source component for the AGN, since the central engine 
is unobscured along the line of sight. There is an emission 
line region west of the cluster which was masked during 
the fits. With no constraints on any of the fit parameters, 
we find that the Sersic component has the largest effec¬ 
tive radius in the F814W band ( R e = 4.56 pc, n = 1.42), 
and the smallest effective radius in F275W ( R e = 1.80 
pc, n = 3.01). We found no obvious trends in the effec¬ 
tive radius or the Sersic index with wavelength. In order 
to examine the degeneracy between the magnitude of the 
Sersic component and the point source component in our 
fits, we repeated the fits in all seven bands with the ef¬ 
fective radius of the Sersic component fixed at 1.80 pc 
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Fig. 17.— Luminosities, effective radii, Sersic indices and axis ratios of the individual Sersic components as a function of filter central 
wavelength for NGC 4244. The cluster consists of an extended, highly flattened, blue component and a compact, spheroidal, red component. 
These components are denoted with a blue short-dashed line and a red solid line, respectively. 


and the Sersic index fixed at 3.01, and then again with 
with the effective radius fixed at 4.56 pc and the Sersic 
index fixed at 1.42. The range of magnitudes obtained 
from these fits gives a rough estimate of the photometric 
uncertainties due to the degeneracy between the profiles 
of the point source component and the Sersic component 
in the fits. Figure [18] shows the SED of both the Sersic 
component and point source component obtained from 
these fits. From single Sersic component fits (with no 
PSF component), we find that the axis ratio monotoni- 
cally increases from 0.55lo‘(n ‘ ln F275W to 0.89^0 oi 1n 
F 153M. _ 

iFilippenko et a! ] (TmMi measured a UV continuum 
slope of = —2 by fitting a power-law model of the 
form f v oc to a HST /Faint Object Spectrograph 
(FOS) spectrum of the AGN over the range 2200 to 
3300 A. We fit a single power-law model to the SED of 
the AGN from the GALFIT decomposition over all seven 
photometric bands using the IDL-based least-squares fit¬ 
ting code MPFIT (jMarkwardtl l2009f ). and measured a 
continuum slope a = —0.8 ±0.2. We also indepen¬ 
dently fit power-law models to the SED in the UV and 
from the optical to the IR, and measured a slope of 
— 1.5 ± 0.6 between F275W and F 336W, which is con¬ 
sistent with IFilippenko et al.l (jl993[ ) , and a measured a 


slope of —0.3 ± 0.4 between F438W and F153M. How¬ 
ever, our observations only represent a single temporal 
snapshot of the state of the AGN, and past observations 
have shown that the optical continuum sha pe of t he AGN 
sometimes undergoes rapid changes. [Lira et all ( 1999 ) 
found that the optical continuum slope changed from 
—2 to 0 over a period of just six months between July 
1996 and January 1997. The SED of the AGN from the 
GALFIT decomposition has an irregular shape, and the 
slope between the A?, V, and I bands is not smooth, sug¬ 
gesting that the fits are strongly affected by degeneracy 
between the AGN and compact cluster component. The 
photometric uncertainties from the GALFIT decomposi¬ 
tion are quite large (^ 0.4 mag for point source compo¬ 
nent and ~ 0.2 mag for the Sersic component in F814W), 
and improving this decomposition will require imaging of 
the NC with higher spatial resolution than what is cur¬ 
rently available. 

5.10. NGC 7793 

At a distance of 3.4 Mpc (|Pietrzvriski et al1l2010f h the 
nearly face-on galaxy NGC 7793 i s one of t he bri ghtest 
galaxies in the Sculptor Group. IWalcher et al.l ( 2005 ) 
measured a stellar velocity dispersion of 24.6 km s -1 
and a mass of 10 6,89 M e for the NC. NGC 7793 was 
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Fig. 18. — SEDs of the individual fit components for NGC 4395. 
A point source component (shown in blue) was used to model the 
AGN and a Sersic component (shown in red) was used to model 
the NC. The solid lines correspond to the SEDs resulting from 
fits performed with no constraints. In addition to the fit with no 
constraints, fits were performed with the effective radius fixed at 
4.56 pc and the Sersic index fixed at 1.42, and with the effective 
radius fixed at 1.80 pc and the Sersic index fixed at 3.01. The 
shaded regions are bounded by the SEDs resulting from taking the 
minimum and maximum magnitude of all three fits in each band. 

included in the iBdker jt ap d2002h #ST/WFPC2 NC 
imaging survey, where an /-band absolute magnitude of 
— 13.64 =b 0.03 mag was measured. Correcting for the 
slightly different distances used, this corresponds to a 
WFC3 F814W absolute magnitude of —13.61 ±0.03 mag, 
which is in very good agreement with our measurement 
of — 13.59lo'o5 ma g- We fit a flat sky to the host galaxy 
and a single Sersic component to the cluster in each 
band. The cluster has a larger extent in the UV than the 
IR, with an effective radius that monotonically decreases 
from 12.45^^83 P c F275W to 7.94 ^i'.o 2 pc in F153M. 
The color gradient is obvious from inspection of the color 
composite image of the cluster shown in Figure 02 The 
Sersic index monotonically increases from 1.39^o‘27 
F275W to 3.29^0 5 g in F153M. Although a single Sersic 
component gives a reasonable fit in all bands, a closer 
inspection of the images reveals a ring-like gap in the 
stellar density, which is particularly prominent in the UV 
bands. This structure can also be seen in the residuals 
of the F814W fit (Figure [5]). 

6 . CONCLUSIONS 

We have analyzed the structure of a sample of ten nu¬ 
clear star clusters in bulgeless, late-type spiral galaxies 
by fitting PSF convolved 2D surface brightness models to 
HST /WFC3 images of the clusters in seven wavebands 
ranging from the near-UV to the near-IR. We find that 
most of the NCs are described remarkably well by a sin¬ 
gle Sersic component, although some show clear evidence 
for multi-component substructure. NGC 2976 has a NC 
which consists of a compact, blue structure, which may 
be a young star cluster and an extended, red population. 


The NC in NGC 4244 contains a red, compact, spheroidal 
component and a blue, extended, flattened component. 
NGC 4395 hosts an unobscured AGN in its NC, which 
was modeled using a point source component. The main 
results of the surface brightness profile fits are as follows: 

1. The clusters exhibit a wide range of structural prop¬ 
erties. In the F814W filter, absolute magnitudes range 
from —11.20 mag to —15.02 mag, effective radii range 
from 1.38 to 8.28 pc, Sersic indices range from 1.63 to 
9.73, and axis ratios range from 0.57 to 0.94. We find 
a much wider range of Sersic indices than any previous 
study, which have generally measured n ~ 1-3. The 
highest Sersic indices we measure (above 6) serve as 
useful indicators of the shape of the wings of the NC 
profile, but do not provide information about the core 
concentration due to the effects of PSF blurring. 

2. We find evidence for spatially segregated stellar 
populations within some clusters. There are color gra¬ 
dients in six of the ten NCs in our sample, indicated by 
a dependence of the measured effective radius on wave¬ 
length. In M33, NGC 247, NGC 300 and NGC 4395, 
we find an increase in the effective radius with wave- 
lengt h. This is cons i stent with general results found 
by iGeorgiev fe Boke~rl ( 2014 ) for a large HST /WFPC2 
archival sample, and indicates the presence of a younger 
population which is more concentrated than the bulk of 
the NC stars. However, we find a general decrease in 
effective radius with wavelength in NGC 3621 and NGC 
7793, which may indicate extended, circumnuclear star 
formation. For NGC 3621, the larger extent of the NC 
in UV bands may be due to scattered UV photons origi¬ 
nating from the obscured Type 2 AGN. 

3. We also find a general trend of increasing roundness 
of the NCs at longer wavelengths, indicated by larger 
best-fit axis ratios in redder bands. There is also a cor¬ 
relation between the axis ratios of the NCs and their 
host galaxies. We conclude that blue disk structures 
aligned with the host galaxy plane are a common fea¬ 
ture of NCs in late-type galaxies, but are more easily 
detected in galaxies that are close to edge-on. 

4. A comparison of the measured colors of the NCs 
with SSP evolutionary tracks clearly shows that in gen¬ 
eral, a mix of young and old SSPs is needed to de¬ 
scribe their stellar populations, consi ste nt wit h pre¬ 
vious spectroscopic studies of NCs (Ro ssa et al.l 12006 : 
iWalcher et, al.ll2006h . The mix of ages is only seen when 
considering colors which span the UV to the IR, and is 
not apparent from optical colors alone. The wide wave¬ 
length coverage of our data provides sensitivity to stellar 
populations with a mix of ages. 

The surface brightness profiles presented in this paper 
will be used in future stellar population modeling and 
dynamical studies of the clusters. Recent studies have 
shown that degeneracies between age, metallicity and ex¬ 
tinction can be broken by fitting SSP models to photo¬ 
metric SEDs of star clusters w ith UV to IR broadband 
coverage (jde Meulenaer et al.112014 ). Our broad wave¬ 
length coverage will allow us to place tighter constraints 
on the star formation histories of the NCs than previ¬ 
ous studies. By performing pixel by pixel SED fits to 
the surface brightness models presented in this paper, 
we can further search for population gradients within 
the NCs and constrain their stellar mass profiles. Stellar 
mass profiles, along with IFU stellar kinematic data of 
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the clusters will be key ingredients in future dynamical 
studies of NCs, which may uncover a previously unseen 
population of intermediate-mass black holes within the 
clusters. 
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